Design method of dynamical decoupling sequences integrated with optimal control 

theory 
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A method for synthesizing dynamical decoupling (DD) sequences is presented, which can tailor 
these sequences to a given set of qubits, environments, instruments, and available resources using 
partial information of the system. The key concept behind the generation of the DD sequences 
involves not only extricating the strong dependence on the coupling strengths according to the 
"optimal control," but also exploiting the "refocus" technique used conventionally to obtain DD 
sequences. The concept is a generalized one that integrates optimal control and designing of DD 
sequences. 

PACS numbers: 03.67.Pp, 03.65.Yz, 82.56.Jn 



I. INTRODUCTION 

In quantum information processing (QIP), dynamical 
decoupling (DD) [J i] is a useful tool for suppressing 
the decoherence resulting from multiple couplings among 
qubits and their environments. Potential applications 
of dynamical decoupling to quantum memory |3| , quan- 
tum computation (QC) [3-0] , and fault-tolerant quan- 
tum computation [8| have been proposed recently. Trains 
of successive pulses such as concatenated DD, Uhrig 
DD (UDD), and quadratic DD (QDD) [1 ©, [ld| have 
been proposed and demonstrated successfully to elimi- 
nate unwanted couplings (Il] - fl3j without knowledge of 
surrounding environment, e.g., its coupling strengths or 
structures. These methods are useful for experiments in 
which system details are largely unavailable. In prac- 
tical experiments, however, theoretical approximations 
such as the delta pulse approximation result in the fi- 
nite pulse width problem [14j . and the requirement of 
an infinite frequency bandwidth for a rectangular pulse 
causes deformation at the leading and trailing edges of 
the pulse [la, LL6|. Both these practical resource limi- 
tations cause serious bottlenecks when using pulse DD 
sequences. In contrast, numerically optimized pulse se- 
quences [lTj-Qjl make use of the complete qubit- and 
environment-related information. However, these meth- 
ods are attractive [2(| only when one can obtain details of 
a qubit-bath system, although spectroscopy and genera- 
tion of pulse sequences become difficult as the system size 
increases. Hence, when using partial system information, 
we can expect better performance from DD sequences. 
Previous studies that have used this approach include 
refs. [U [2~3 ], in which the spectra of the surrounding 
environment were reflected on their energy-constrained 
decoupling waveforms. 

In this paper, we propose a method to numerically 
synthesize DD sequences under realistic resource limita- 
tions such that they can be tailored to a system based 



on available knowledge about qubits, environments, and 
instruments. In the two examples given in this paper, we 
generate bandwidth- and energy-constrained DD wave- 
forms that are free from the finite pulse width problem 
and are robust to waveform deformation arising from in- 
struments [16]. The paper is organized as follows: In 
Sec. UU we describe the basic concept of our DD design. 
In Sec. IIII1 we present a pragmatic example in which a 
single qubit is coupled to many two-level systems, and 
where the control field has imperfections due to instru- 
ments. In Sec. IIV1 we apply our method to a multiqubit 
system where the qubits couple to each other. In Sec. [V] 
we discuss the performance of DD sequences generated 
under resource restrictions in terms of energy and band- 
width. The appendices describe how the average Hamil- 
tonians are calculated using Floquet Hamiltonian theory 
in the examples. 



II. BASIC CONCEPT 

Let "Ho be the Hamiltonian of a system under consider- 
ation, and V(t) be a time-dependent external perturba- 
tion, represented in an appropriate rotating frame. We 
use H = 1 units throughout. The system Hamiltonian Hq 
can be expressed as 



(1) 
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using a set of orthonormal Hermitian operators 
{0 Q |0 Q = 0t, Tr[8 Q 9,g] = S a p} and their coefficients 
{h a }. For example, a system Hamiltonian that de- 
scribes the coupling between a single qubits and the sur- 
rounding environment can be expanded with operators 
9 Q = \o a ® B a (a — x,y,z), where a a are the qubit's 
Pauli operators and B a are bath operators. Consider Hq 
in the interaction frame of V(t) using unitary operator 
Uv(t) = Texp(i f Q V(t')dt'), where i = \f—l and T is 
the time-ordering operator. The Hamiltonian of Eq. (fTJ) 
in the frame of V{t) then becomes the modulated Hamil- 
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tonian Ho(t) as follows: 



H (t) = Uv(t)HoUl(t) = EM E^W^ 9 ' ( 2 ) 



where tjy (t)@ a Uy(t) = J2p c /3a{t)Qp- In this process, 
the effect of finite pulse width is naturally incorporated 
into Cff a (t). Delta pulse approximation of V(t) can be 



also treated as discontinuities of cp a (i). We call cp a d) a 
system- modulation matrix, as used in ref. (22[. 

This system can be stroboscopically "time-suspended," 
that is, decoupled from environment, if there is an exter- 
nal perturbation V(t) such that the average Hamiltonian 

H cS of H (t) vanishes @, [TJ, [Hj], and if the period T 

of H e s satisfies ||7io?1| < 1 so that Magnus expansion 

of the evolution operator converges. Suppose Hod) has 
T-periodicity; then, the average Hamiltonians @, [23| up 
to the first order become 



^(o) i 



H cS = - I Hod')dt' = J2 h <* ( f L c ^ dt ' I 9 ^ 



T 

— (!) 1 \ - 

^eff = 2iT haha ' 1^ 

a, a' P<P' 



a, P 



dt 2 / dti {Cp> a >(t2)cp a dl) - C0> a dl)Cf3 a >d2)} [0/3', 8/9]- 
JO I 



(3) 
(4) 



We can obtain DD sequences Vd) by minimizing the co- 
efficients of the operators Qp and [©#', s], i.e., the terms 
within the large curly brackets in Eqs. © and (Q}. Given 
that the external perturbation Vd) does not always guar- 
antee the T-periodicity of "Ho(i), the system- modulation 
matrix cp a d) should be constrained such that cp a d) are 
continuous at t = 0, T, • • •. 

For numerical minimization, Vd) should be parame- 
terized with pulse timings, flip angles, phases, and so on. 
When the parameters are denoted by £ = {(i,C2, • • •}, 
the system-modulation matrix is also parameterized such 
that Cp a d) = cp a d] 0- A cos t function for the numerical 
minimization is constructed such that 



(i) 

a'a,P'P 



(0 



(5) 



a,p 



a'a,p'<P 



together with a penalty function of the continuity of 

cp a d) at t = 0, T, 



$p(0 = El c /9«( T ;0-c^(o;c)i 2 

a,P 



where 



cp a d';t)dt' 







(0) 



(7) 



®Xp>p(0 = t 



T rt 2 

dt2 / dti 
o Jo 

X {c / 9/ Q .'(t 2 ;C)e/9a(*i;C) 
-Cp> a dl'X)cp a -d2\ C)} 



(8) 



Consideration of higher-order perturbations would re- 
sult in better decoupling performance; however, doing so 



would involve cumbersome hand calculations for deriving 
average Hamiltonians. 

The synthesis of DD sequences is then reduced to a 
nonlinear minimization problem: 



mjn$(C) subjected to $p(C) = 0, 



(9) 



which can be numerically performed on a modern digital 
computer. 

This method is basically the same as that given by 
Eq. (fO) in [f8| except that the cost function is con- 
structed in Lie algebra. However, this difference allows us 
to choose among sequences of DD type, optimal control 
type [17h1j1> or a combination of both these types. In 
general, the expansion coefficients {h a } in Eq. ([TJ con- 
tain the details of the system, e.g., coupling strengths 
and frequency shifts, which can be obtained through 
spectroscopy. In the synthesis of DD type sequences, 
which do not require knowledge of the qubit's surround- 
ings, {h a } should be excluded from the cost functions in 
Eqs. §7$ and §5$ so that the generated sequences do not 
depend on them. For the case where {h a } are known, se- 
quences can be synthesized with their help. If we explic- 
itly include {h a } in a cost function, for example, Eq. ([7]) 
becomes 



^ 0) (C)=T- 1 



E ft « / cpad'-xw 

Jo 



(10) 



and the obtained sequences suppress decoherence using 
that information. In the limit of ||"HT|| — > 0, or in the 
limit of no — > 0, where uq is the maximum order of an 
average Hamiltonian, this method corresponds to uni- 
tary matrix based optimal control. When {h a } are par- 
tially known, they can be incorporated into the cost func- 
tion to produce a hybrid type of DD and optimal con- 
trol. This would be advantageous when we know the 
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qubit-qubit interaction strengths that are to be used for 
gate operations, but do not know the strengths of qubit- 
environment couplings that cause decoherence. 

Furthermore, cost functions in Lie algebra alleviate the 
growth of computational complexity as a system becomes 
larger. Given that most operators in Lie algebra, which 
stand for many-body interactions, barely appear, there 
is a dramatic decrease in computational costs associated 
with numerical optimization for a system Hamiltonian 
that is simple and highly symmetric, i.e., it is described 
with a small number of operators in Lie algebra. We 
give a paradigmatic explanation in which a Hamiltonian 
Ho = Y^k=o9^&z,k dictates a system, where a z ^ is a 
A;th qubit's Pauli z-operator, {g^} are inhomogeneously 
broadened energy shifts, and n is the number of qubits. 
In optimal control type, to obtain a DD sequence (time- 
suspension sequence), we have to compute 2™ x 2" ma- 
trices, which are intractable for large n. However, a well- 
known method that flips the sign of a Zj k in the frame 
of V(t) by applying consecutive tt pulses to the qubits, 
such as CPMG [2J|, refocuses the inhomogeneity. Our 
method automates such a flipping procedure for a given 
system. 

Note that while synthesizing DD sequences, Ho should 
be modulated by Uy(t). This requirement is equivalent to 
the bracket generation condition mentioned in ref. (l8j . 

In the following two examples, in order to impose a 
frequency bandwidth limitation on V(t) = ^ Q v a (t)Q a , 
we use Fourier coefficients of 
v 

v a (t) — v aj7l sin(2n7rf/T) (11) 

n=l 

as the parameters of V(t), so that C = {v a ,n}- The Flo- 
quet average Hamiltonian theory [25[ is useful to calcu- 
late the average Hamiltonian when a system-modulation 
matrix is T-periodic and expressed as the Fourier expan- 
sion c/3 a (t) — J2 n c Pa,n exp (2iimt/T). A brief explana- 
tion of this calculation is given in Appendix [A] and the 
cost functions in the following examples are built using 
the theory. 

III. DEPHASING PROBLEM 

Let us consider the qubit dephasing problem. A single 
qubit is coupled to a bath consisting of many two-level 
systems (TLSs) without energy relaxation. When there 
exist inevitable pulse errors caused by instruments, an 
experimenter might need a DD sequence that is robust 
to these errors. In this example, we show how to incor- 
porate Tycko's composite pulses [26[ into DD designing. 
In ref. [261 ] , a Hamiltonian of an external perturbation to 
control the qubits is described as V(t) + where 
V(t) denotes an ideal operation to the qubit and H-E^t) 
represents systematic errors due to the instruments. Self- 
compensating composite pulses are built so that the ideal 
pulsing of V(t) decouples error Hamiltonian 



TABLE I. Parameters of the DD sequence generated in 

Sec. ma 



n 


Vx,nT /it 




1 


-0.7030256 


-3.6201768 


2 


3.3281747 


3.8753985 


•> 
o 


1 1 ^Qflfl77 


1 9^1 1 Q1 Q 


4 


2.9375301 


-0.2998110 


5 


-1.8758792 


3.1170274 


6 


1.7478474 


0.3956137 


7 


5.6966577 


-0.3593987 


8 


-0.5452435 


-3.5266063 


9 


4.0826786 


2.4900307 



Let the Hamiltonian of the system be Ho = Hq-b + 
W.E(t) and the ideal external perturbation be V(t) — 
v x (t)a x + v y{t)&yi where the qubit-bath interaction is 
Wq-b = Y^rk 9 ^z^z,k and the pulse error term is 
He = £\v x (t)a x +£sv y (t)a y +£2Vy(t)a x +EiV x {t)a y . Here, 
a a and a a ,k denote Pauli operators of the qubit and the 
kth TLS, respectively. e\ and £3 denote the error am- 
plitudes of flip angle error, while £2 and £4 denote the 
error amplitudes of phase orthogonality error arising, for 
example, from microwave IQ mixers or hybrid circuits. 
The external perturbation V(t) modulates the Hamilto- 
nian Ho in its interaction frame, so that 

k 

+ y^{£iv x (t)cp x {t) + e 3 v y (t)cp y (t) 

+£ 2 Vy(t)Cp x (t) + EiVxi^Cpy^d-f}, 

where cp a {t)(jp = Uv(t)a a U v {t). In this example, 
the zeroth-order cost function is defined as 





2 


T 2 


/ cp z (t)dt 


+E 


/ v x (t)cp x (t)dt 


Jo 





Jo 



and the first-order cost function is also defined in the 
same manner (See Appendices [5] and [C] for detailed 
derivations). Exclusion of and £j (J = 1, • • • , 4) from 
the overall cost function gives DD type pulse sequences 
with a tolerance to pulse imperfections. DD sequences 
generated with this method can be tailored to the sys- 
tem by multiplying deliberate weights to terms in the 
cost function. We minimized the cost function by op- 
timizing £ = {!)„.„} using a steepest descent method 
in combination with a genetic algorithm. The gener- 
ated waveform expressed in Fourier series form, v a (t) — 
Xm=i Va > n sin(27rni/T), as listed in Table is shown in 
the top panel of Fig. Q] 

We evaluated the performance of our synthesized DD 
sequence, UDD12, and QDD3 against the flip angle error 
A/3 = £1 = £3 and the phase orthogonality error Aip 
= t&x\~ l (£2V y {t)/v x (t)) on a single qubit system coupled 
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FIG. 1. (color online). Decoupling fidelities against (a) flip 
angle error and (b) phase orthogonality error for a single qubit 
coupled to four bath TLSs with UDD12 and QDD3 sequences, 
and the sequence generated in this example. The waveforms 
of these sequences are displayed in the top panel of the figure. 



to four bath TLSs, with the fidelity of quantum gate C 
defined as 

F(C) = min w || VWMy/Cmmitr. (12) 

The number of pulses and the pulse widths in UDD 
and QDD sequences were chosen such that the total 
amounts of the applied energy and the peak ampli- 
tudes were nearly equal to those of the generated wave- 
form. If the energy of a DD sequence is defined as 



Jo J2a=x. y V a( t ) tne energies per cycle of the gen- 
erated DD sequence, UDD12, and QDD3 amount to 
11.50 vr/T, 12 tt/T, and 16 tt/T, respectively. The 
peak amplitudes max tg r T ] w 2 (t) T/2ir of the DD 

sequences were limited to 10, 20, and 22 for the syn- 
thesized DD waveform, UDD12, and QDD3, respectively. 
Given that there is minimal interpulse spacing in UDD 
and QDD, their peak amplitudes were not equal to that of 
the synthesized DD sequence. The coupling strengths be- 
tween the qubit and TLSs were randomly generated; we 
used {g {k) T/Tr} 4 k=1 = { 0.0338264, -0.0906347, 0.0014495, 

0.0740895 }, so that ||# T||/27r ~ 0.1. 

As the flip angle errors are accumulated at each tt 
pulse, UDD is sensitive to the flip angle error (Fig.[TJa)); 
on the other hand, whereas QDD cancels the flip angle 
error by inserting a 90° phase-shifted tt pulse. However, 
as seen from Fig. [IJb) , QDD cannot compensate for the 
phase orthogonality error. Figures Ufa) and (b) show 
that the generated sequence is robust against the pulse 
errors caused by the instruments. This example shows 
that a DD sequence can be optimized to be robust against 
known systematic errors, which are small but accumulate 
after millions of cycles in actual experiments. A variation 
of this example, where there is crosstalk between control 
pulses, can be managed in the same manner. 

Another variation is that the behavior of the surround- 
ing bath TLSs is known, the Hamiltonian of which is de- 
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FIG. 2. (color online). DD performance as evaluated for n- 
qubit chain system coupled through dipolar interactions with 
(a) optimized pulse sequence and (b) MREV16 sequence. Fi- 
delity F is normalized so that F = 1 — (1 — F)/\\Hq-q \\ 
because H'Hq-qH becomes larger as the number of qubits in- 
creases. 



noted by Hb, e.g., the principal axes of TLSs are tilted 
from the z axis so that [Hq-b^b] 0- I n this case, 
taking the bath Hamiltonian of the TLSs into account, 
we can obtain a DD sequence that overcomes degrada- 
tion due to H.b by minimizing higher-order cost functions 
$( n \ at least up to the second order. 



IV. MULTIQUBIT COUPLING 

Consider a system that consists of a finite number 
of qubits coupled to each other. Access to the qubits 
is assumed to be restricted to collective manipulation, 
the external perturbation of which is given by V(t) = 
v x(t) J2k °x,k + v y {t) J2k &y,k- Each qubit is considered 
as being connected through a qubit-qubit Hamiltonian, 

n = n Q -Q = E dik '' k) 

k'<k 

X (PzM'GzM — Vx.k'Ox.k/Z - Oy.k'Oy.k/'Z) . (13) 

In the interaction frame of V(t), we obtain a modulated 
Hamiltonian, 



(14) 



where 77/37 (i) = cp z (t)c 7Z (t) - cp x (t)c lx {t) /2 - 
cp y (t)c ly (t) /2. We defined a cost function for the syn- 
thesis of v x (t) and v y (t) in a manner similar to that in 
Sec. IIII1 excluding d^ k ^ from the cost function to man- 
age the arbitrary strengths of the interactions and con- 
sidering ||Hq_qT|| < 1. The synthesized waveforms ex- 
pressed in a Fourier series are listed in Table [TT1 

We evaluated the synthesized DD waveform by simu- 
lating a one-dimensional qubit chain by considering cou- 
plings up to the second nearest qubits, and compared its 
performance with that of MREV16 [27J, which is known 
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TABLE II. Parameters of the DD sequence in Sec. IIVI 



n 


Vx,nT/-K 


Vy,„T/lY 


1 


-4.8892576 


-2.6291726 


2 


-3.1490576 


-3.4112889 


3 


-14.317448 


-1.7326439 


4 


-0.0929321 


4.2805093 


5 


6.8394959 


-3.7925374 


6 


-0.6645375 


-2.3678092 


7 


0.3344480 


-2.5797746 


8 


-1.5042059 


-1.9232075 


9 


2.3863574 


-4.2795712 



to be powerful for this type of coupling. The Hamilto- 
nian of the qubit chain was denned using the coupling 
strengths in Eq. (fl"5|) . 

{ir/T (k' = k + l) 
n/8T (k' = k + 2) . (15) 
(otherwise) 

Because it was difficult for us to calculate the gate fi- 
delity defined in Eq. (fT2|) for several qubits, we used the 
state fidelity for typical quantum states and the trace fi- 
delity defined as F(C) = Tr[V C^C] instead, where C is 
the unitary representation of C. Four states were used 
and defined as follows: The coherent spin state directed 
toward the x-axis is the eigenstate of the collective x- 
operator a x ^ that has the maximum eigenvalue. The 
GHZ state is |GHZ) = |0)® n + |1)®", where n is the 
number of qubits. The maximally entangled state is 

|MES) = 2~ n/2 J2Zo I*)- Thc Dicke state is one of thc 
eigenstates of the collective z-operator and the total spin 
operator a ■ a that have the a z eigenvalue closest to zero 
and the maximum a ■ a eigenvalue. The energy and peak 
amplitude of this DD sequence were restricted to 13.8 
tt/T and 20 tt/T so that they were almost equal to those 
of MREV16 (16 tt/T and 20 tt/T, respectively). 

Figure [5] shows the fidelity per cycle against the num- 
ber of the qubits in the system. The synthesized se- 
quence effectively decouples qubit-qubit interactions for 
any state available in this example, whereas MREV16 
shows considerable degradation in the nonclassical states 
owing to cooperative destruction caused by the finite 
width of the pulses and the qubit-qubit interactions. 
Given that the numerically synthesized sequence is op- 
timized to isotropically suppress any operator coefficient 
of up to the first order, experimentally, it works better 
for the nonclassical states, which are necessary for QIP 
experiments. 



V. PERFORMANCE VERSUS RESOURCES 

Our examples, with smoothly modulated wave DD se- 
quences, are free from the finite pulse width problem 




FIG. 3. (color online). Performance under resource restric- 
tions in terms of amplitude of ||V r (f)T|| and bandwidth of 
V(i). Regions in which the maximum amplitude or the band- 
width is excessively supplied are labeled as I and II, respec- 
tively, and the region where both resources are sufficiently 
supplied is labeled as III. All axes are in log scale. 



and are robust against distortion in the waveforms under 
the finite frequency bandwidth limitation. Available re- 
sources can be specified in terms of the maximum control 
amplitude || V(i)T|| and frequency bandwidth of V(t), so 
that DD sequence synthesis is performed within the avail- 
able resources. We found that there is a close relationship 
among the DD performance, available bandwidth, and 
peak amplitude. Figure [3] shows the performance of DD 
sequences used in Sec. IIVI under the resource restrictions. 
The root mean square of the operator coefficients in an 
average Hamiltonian up to the first order v <J>(°) + <J>W 
was used as a measure of performance because this value 
indicates the suppression ratio of a coupling strength at 
H^o^ll ~ 1. We imposed frequency bandwidth limitation 
on V(t) using p in Eq. (fPT|) . Bandwidth of waveforms 
/bw was defined as 2p/T. 

The figure clearly shows that the performance of the 
generated DD sequences monotonically improves with 
an increase in the available resources; however, it also 
suggests that there is a trade-off between the perfor- 
mance and available resources. As seen from thc fig ure, 
converged values of the cost function v $(°) + $W are 
closely related to the available resources. This can also 
be seen in ref . [22j , where the performance of synthesized 
DDs depends on available amplitude of V(t). Although 
the converged values are not always the global minima of 
<5(C) (because of the existence of local minima in nonlin- 
ear minimization), the figure suggests that available re- 
sources bound the performance of dynamical decoupling 
sequences. 

The figure shows another relationship between the re- 
sources: the decoupling efficiency drops when either of 
the resources is oversupplied. Region I in Fig. [3] is as- 
cribed to over driving, which denotes lack of bandwidth 
of V(t) to utilize the large amount of driving energy 
for controlling the qubits. This behavior can be seen 
in Fig. 2] (I), where a part of a system- modulation ma- 
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FIG. 4. (color online). System-modulation matrices c ax (t) (a 
= x, y, z) in Regions I, II, and III, represented in Cartesian 
coordinates. c ax {i) are displayed as functions of t £ [0, T]. 

trix c ax (t) (a = x, y, z) in the region is represented as a 
trajectory in Cartesian coordinates. Given that the tra- 
jectory needs to have loops to consume the excess driv- 
ing energy, the trajectory of the system-modulation ma- 
trix is constrained from moving efficiently to decouple, 
and consequently, the performance drops. Region II in 
Fig.[3]is ascribed to over modulation of V(t), wherein the 
amount of energy supplied for driving the qubits is insuf- 
ficient (Fig. @] (II)). When both resources are sufficiently 
supplied to drive and modulate, as in the region as III in 
Fig. [3l the system-modulation matrix in the region effi- 
ciently averages out unwanted interactions (Fig. 0] (III)). 



VI. CONCLUSION 

We presented a general framework to numerically syn- 
thesize DD sequences, which can tailor them to a given 
set of qubits, environment, and realistic resource limi- 
tations. To this end, a cost function for the numerical 
synthesis of the DD sequences was built using system- 
modulation matrices of operators in Lie algebra, instead 
of a unitary propagator, in order to fully exploit the sym- 
metries that exist in qubits and environments. The ef- 
fects of finite pulse width are naturally incorporated in 
the system-modulation matrix. We presented two exam- 
ples to demonstrate the robustness of the generated DD 
waveform to flip angle and phase orthogonality errors, re- 
gardless of the error amplitude, and the capability of the 
DD waveform eliminating unwanted qubit-qubit interac- 
tions, regardless of the number of qubits with protected 
nonclassical qubit states. The trade-off between perfor- 
mance and resource limitations was numerically stud- 
ied. In the near future, we intend to rearrange a system 
Hamiltonian with the help of an external perturbation 
as introduced in [ijj], which is also applicable to QC and 
QIP experiments. 
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and a a , cr a ,k (a = x,y, z) denote Pauli operators of the 
qubit and the fcth TLS, respectively. A unitary operator 

Uy{t) — Texp(j J„ V(t')dt') modulates the Hamiltonian 
Ho in the interaction frame of V(t) such that 



,k e 



2-nint/T 



(B3) 



T~iE{t) — < £i Vx,n—ri c ax,n' 
a,n v n' 

H - &2 ^ ^ ^y,n— n'Cax,n' 

n' 

~t~ £3 ^ ^ ^n—n' Cay ,n' 



£4 



x,n—n'*-'ay,n' f u a 



a a e 



2-jrint/T 



(B4) 



Appendix A: Floquet Hamiltonian Theory 

A periodic time-dependent Hamiltonian T-L{t) can be 
averaged over a period T using an approach based on 
the Floquet Theory. Consider a Fourier expansion of 

H(t) such that 



H(t)=J2 #F ex P (2wrnt/T) . 



(Al) 



As stated in [25|, average Hamiltonians up to second- 
order are given by 



^(0) 

^eff 
^(1) 

?4ff 



^(2) 



n 



(0) 
F > 



(A2) 
(A3) 



T 2 

— V V 

127T 2 ^ ^ n(n + n') 



1 + (1/2)^,0 



(A4) 



Appendix B: Derivation of average Hamiltonian for 
dephasing problem 

The Hamiltonian of the system Hq is considered as 
Hq = H.q-b + He in the main text, where 

Kq-b = ^2g (k) v z v z ,k, (Bi) 

k 

He = e\v x (t)a x + e 3 v y (t)a y 

+ £2V y (t)a x + 6iV x (t)a y , (B2) 



where v a ^ n = T 1 J Q v a (t) exp(— 2iTmt/T) and cp a (t) — 
Y,n c Pa,nexp(2nint/T) = Tr[U v a a Ulap]/Ti[&^]. If we 
write 



-7V ^ax.rv • 



u a,n ~ / j u x,n-n'^ax,n' , u ' a n — v y,n- 
n' n' 

^L,n = ^ ] Vy,n-n'Cgy,n' , <^L,n = ^ ] v x,n—n' c ay,n' 



the Fourier series of the modulated Hamiltonians are 



Hq_ b y — ^ I 9^ ^ Caz,n&z&z,k: 



nj{n) _^ ii) ~ 



(B5) 
(B6) 



The zeroth-order average Hamiltonian is easily obtained 



.(0) 



from H cS = T-L^q_ b f + Hg' F> and is given by 

Wcs = 9 ik) c a z,aa z a z , k + £ * d *fi (B7) 



(o) 



a,k a,i 

The first-order average Hamiltonian is obtained from 
Eq. (|A3|) , and is given by 

T7 (1) _ T ST 1 ra/(») ij(~nh 

— 27T ^ 71 Q~ B ' F n E,F> n Q-B,F n E,F J 
n>0 

2tt ^ 7i 1 Q- S > F ' 0- 



E /r7v( n ^ ■}/ 1 7 / 1 " > 17' ■" ' 
|L n -Q-S,F' "-Q-B.fJ "T L n '£;,F' rl -E,F 

n>0 

L n 'Q — B,F' n 'E,F J l^Q-B.F' n -E,FJj 



where 
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Q — B,F' rL Q — B,F 



^a' z.n^az, — n 



<J a > i<J a \<Jz,k' Oz,k 



-ij( n ) -V/(-«) 
n E,F' n E,F 



nj( n ) o>(-™) 
n Q — B,Fi n E.F 



= E^V fe) E c ° 

k' ,k a' ,a 

^ ^ Q ^ 9^ ^ ^ ^ {Ca' z,nCa.z, — n Caz ,n^a.' z , — n) [^"q ; j ^'aj^'z,^' & z.k 
k' ,k a' <a 

= E^^ E ( C "'^« C L,™ - c az ,„c*, 2in ) [dv , a a ]vz,k> v z ,k 

k' ,Aj a'<a 

= ^2g {k } 9 ( - k) Y 2e Q ' Q/3 Im[c QX ,„c*, 2i „] <x z , fe ' a z . k , 

k' ,Aj a ' < a 

a' <a i' < i 

= E E # (fc)£l E 2e a , Q/3 Im[ C ;, Z:n dW, + CaZ)W d«*J ^ 



and e a /3 7 is the Levi-Civita symbol. X) Q '< a represents we used the identity c a p.- n = c* ^ n because c a p(t) is 
the summation of (a', a) = {(a;, y), (x, z), (y, z)}. Here, real. The first-order average Hamiltonian then becomes 



— (!) T i II \ 

HeB = — Y J 9 Y £a ' a f } [ Y - lm {Caz,nC* a , z ^] \ &0 O z M & z . k 
k' ,k a' <a \n>0 / 

+ T - Y e ^ fe ^m^,„«) 

a' <q i \n>0 / 

+ f E E ^ (y ^m<$,„ <C + ^ 

a'<Qi'<2 \n>0 / 

+ v E E 5 (fc) ^ E ( E ~ Im K<s,„ + c QZ ,„ d^; n ] ) ^ 

fe i a'<a \n>0 / 



(B8) 



Appendix C: Derivation of the cost function for 
dephasing problem 

The zeroth-order cost function is defined from 
Eq. (|B7l) as follows: 



* (o) = e {i^oi 2 +u^|e| 2 } 



(Cl) 



a' <a i n>0 
a.'<a>i'<i n>0 
+ ^ E E E^^'^^+^ft 



n>0 



The first order cost function $W is defined such that 
terms within the curly bracket in Eq. ()B8|) become zero, 
and is given by 



(C2) 



a' <a 



E 



c az,nC a r z n \ 



n>0 



where w is the deliberate weight for the pulse errors set 
to 1/100. 



